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To  transfer  energy  from  one  rectangular  wp.veguide  into  another,  the 
two  guides  are  made  parallel  to  each  other,  as  shown  in  Fig.  1-b.  In  the 
case  of  identical  waveguides,  the  coupling  between  them  depends  on  the  dis- 
tance separating  them  and  the  distance  they  have  in  common  along  their  axis. 
To  control  the  transfer  of  power  from  one  guide  to  the  other,  we  have  to 


study  the  effect  of  these  two  parameters. 

Miller  [2]  studied  the  coupling  between  transmission  lines.  Marcatili 
adopted  his  solution  and  found  the  coupling  coefficients  for  the  two  identical 
rectangular  waveguides  [3] . However,  Miller  used  approximations  which  are 
valid  only  when  the  separation  distance  between  the  two  guides  is  extremely 
large  and,  consequently,  the  coupling  is  very  weak. 

The  purpose  of  our  study  is  to  develop  a general  theory  for  the  coupling 


between  two  identical  rectangular  waveguides.  We  will  present  a different 


and  more  general  method  than  that  of  Miller  and  Marcatili,  that  is  valid  for 
a strong  coupling,  i.e.  for  an  arbitrary  separation  distance.  In  our  cal- 
culation, we  consider  the  reflection  and  transmission  coefficients  between 
the  two  guides  and  use  their  exact  values  in  our  differential  equations  for 
the  power  transfer  between  the  two  guides.  We  prove  that  these  coefficients 
are  complex  quantities  with  a nonzero  real  part.  Unlike  us,  Miller  defined 
a coupling  coefficient  which,  in  his  solution,  turns  out  to  be  a purely 
imaginary  quantity. 

The  geometrical  description  of  the  system  of  two  guides  and  the  state- 
ment of  our  problem  are  given  in  Chapter  I. 

In  Chapter  II,  we  summarize  Marcatili's  study  [1]  of  the  single  rectan- 
gular waveguide.  However,  in  our  presentation  we  use  the  concept  of  ray 
optics;  since  the  electric  and  magnetic  fields  are  treated  as  incident, 
reflected  and  transmitted  waves,  we  use  exponential  forms  for  the  fields 
inside  the  waveguide  instead  of  the  sinusoidal  forms  used  by  Marcatili. 

In  Chapter  III,  we  study  the  system  of  two  waveguides.  In  the  coupling 
region  (0  - z - L),  we  calculate  the  reflection  and  transmission  coefficients 
at  the  inner  boundaries  and  find  the  eigenvalue  equations  for  modes.  The 
Miller-Marcatili  coupling  solution  is  also  discussed  in  this  chapter. 

Chapter  IV  is  reserved  for  a presentation  of  our  method  for  coupling 
calculations  between  the  two  guides,  with  a discussion  and  comparison  with 
the  Miller-Marcatili  method. 

A presentation  of  the  main  results  of  our  analysis  was  given  in  [4] . 


CHAPTER  I 


PROBLEM  PRESENTATION  AND  GEOMETRICAL  DESCRIPTION 

Rectangular  waveguides  are  made  of  dielectric  material  immersed  In  sub- 
strates of  different  dielectric  materials  with  refractive  indexes  smaller 
than  that  of  the  guide.  If  this  condition  is  obeyed,  there  are  angles  of 
incidence  at  the  boundary  core  substrate  for  which  total  internal  reflection 
occurs,  as  a result  of  which  the  light  remains  contained  in  the  waveguide. 
Figs.  1-a  and  1-b  show  the  system  of  single  rectangular  waveguide  and  two 
parallel  rectangular  waveguides,  respectively. 

For  the  system  of  two  waveguides,  the  guides  are  assumed  to  be  of 
identical  cross  section  and  made  of  the  same  material  with  refractive  index 
n^^,  while  the  refractive  index  of  the  substrate  is  n^  < n^ . Also,  the  wave- 
guides are  assumed  to  be  buried  inside  the  substrate  at  a sufficiently  deep 
distance,  so  that  we  can  neglect  the  effect  of  the  boundary  air-substrate. 
Fig.  2 shows  the  two  guides  in  the  transverse  plane  (x,y).  The  two  rectan- 
gular cross  sections  are  surrounded  by  the  materials  of  the  substrate.  We 
consider  the  more  general  case  where  the  materials  on  all  four  sides  of  the 
guides  have  different  refractive  indexes,  as  shown  in  Fig.  2. 

We  define: 

Region  1 - where  the  refractive  index  n = n^. 

Region  2 - where  the  refractive  index  n = n2» 

Region  3 - where  the  refractive  index  n = n^, 

Region  4 - where  the  refractive  index  n = n^,  and 

Region  5 - where  the  refractive  index  n = n^. 

The  separation  distance  between  the  two  guides  is  c,  the  width  of  each  guide 
in  the  y-direction  is  b,  and  the  height  in  the  x-direction  is  a. 


Fig.  3 shows  the  two  guides  in  the  y-z  plane.  In  the  region  (z  5 0) 
we  have  only  one  guide  (we  call  it  guide  I);  in  the  region  (z  i L)  we  have 
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the  other  guide  (we  call  it  guide  II);  in  the  region  (0  $ z ^ L)  both  guides 
are  present  (we  call  it  coupling  region).  In  the  coupling  region,  a transfer 
of  power  occurs  back  and  forth  between  the  two  guides.  We  call  L the 
coupling  distance. 

We  assume  that  we  have  a wave  inside  guide  I traveling  in  the  positive 
z-direction,  and  that  the  condition  of  total  internal  reflection  is  satis- 
fied, so  that  all  the  optical  power  is  contained  in  guide  I. 

When  the  wave  reaches  the  boundar>  z = 0 where  guide  II  begins,  we  have 
a transfer  of  power  from  guide  I to  guide  II.  Traveling  waves  occur  in  each 
of  the  two  guides  with  amplitudes  changing  along  the  z-axis.  At  the  boundary 
z = L,  the  wave  in  guide  I is  partially  reflected,  and  the  wave  in  guide  II 
continues  to  propagate  in  the  positive  z-direction.  Because  of  the  fact 

that  guide  I and  guide  II  are  identical,  the  guided  modes  in  guide  II  (for 

z ^ L)  are  the  same  as  the  guided  modes  in  guide  I (for  z 5 0) . As  a conse- 
quence, the  wave  in  guide  II  (for  z - L)  will  carry  the  original  incident 

mode  and  the  condition  of  total  internal  reflection  will  also  be  satisfied. 

Inside  guide  I,  we  call  incident  wave  the  wave  reaching  the  boundary 
(z  = 0)  from  the  left;  inside  guide  II,  the  wave  traveling  to  the  right  of 
the  boundary  z = L is  called  the  transmitted  wave. 

We  assume  that  we  know  the  power  associated  with  the  incident  wave, 
and  we  want  to  determine  the  transmitted  power  in  function  of  the  sepa- 
ration distance  c and  of  the  coupling  distance  L. 

In  our  study,  we  consider  the  case  where  the  power  associated  with  the 
partially  reflected  wave  at  the  boundary  z = L is  very  small  compared  to  the 
transmitted  power;  hence,  we  ignore  it,  which  means  that  we  will  be  dealing 


only  with  waves  traveling  in  the  positive  z-direction. 


To  determine  the  electric  and  magnetic  fields  and  to  find  the  guided 
modes,  we  solve  Maxwell's  equations  in  three  different  regions  (z  - 0, 

0 - z - L and  z - L) , and  impose  the  boundary  conditions.  Then,  the  trans- 
mitted power  in  function  of  the  incident  power  is  calculated. 

In  the  regions  (z  ^ 0,  z ^ L)  we  have  to  study  the  single  rectangular 
waveguide.  As  we  mentioned  in  the  introduction,  this  study  was  done  by 
Marcatili  [1],  and  will  be  summarized  in  Chapter  II.  The  electric  and 
magnetic  fields  will  be  written  down  in  exponential  form  inside  the  guiding 
core;  thus,  the  form  of  our  equations  will  be  different  from  Marcatili's. 

In  the  region  (0  - z ^ L),  we  have  two  guides,  and  the  solution  will  be 
worked  out  in  Chapters  III  and  IV. 


iFim 
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CHAPTER  II 

RECTANGULAR  DIELECTRIC  WAVEGUIDE 


We  will  study  the  rectangular  waveguide  shown  in  Fig.  3 in  the  region 
(z  - 0).  In  the  x-y  plane,  the  waveguide  has  the  cross  section  represented 
by  the  rectangle  at  the  left  of  the  x-axis  of  Fig.  2. 

We  assume  that  the  time  and  z-dependence  are  given  by  the  factor: 


e-i(6pZ  - o)t) 


(1) 


2 

k = Cl) 


2tt  . 


Where  u is  the  angular  frequency,  it  is  given  by  tne  relation: 

where  k = -^  is  the  wave  vector  in  free  space.  The  quantity  6,,  is  the  propa- 
A 0 

gation  constant  along  the  z-axis  in  the  region  (z  ^ 0) . The  factor  of  eq. 

(1)  will  be  omitted  in  the  following. 

We  will  consider  the  case  where: 

(3) 

This  corresponds  to  modes  far  from  cutoff;  thus,  almost  all  the  power  is 
contained  inside  the  core  of  the  waveguide.  In  this  case  we  can  neglect  the 
fields  in  the  shaded  areas  of  Fig.  2. 

2.1  Guided  Modes  of  the  Rectangular  Waveguide. 

To  get  a complete  discription  of  the  guided  modes  inside  the  waveguide, 
we  must  solve  Maxwell's  equations: 


$0  = n^k. 


VXH  = f . 

= -“o  I- 


(4) 

(5) 


where  the  operator V is  given  by: 


J 


V = ■r—  e + -7—  e + -j—  e , (6) 

6x  X 6y  y 6z  z* 

and  n is  the  refractive  index. 

Considering  relation  (1),  we  can  write  the  transverse  field  components 
in  function  of  the  longitudinal  components  as: 


Ejj  =-  J 


u -Lfe  ^ 1 

Hx-  7TLP*\^  j j 


u ^fPo— J 

- ^2  Lro^^  j J 


Here  j indicates  the  region  number  in  the  x-y  plane,  and: 

‘'j'  = "j'"'  ■ 

The  longitudinal  components  E and  H must  satisfy  the  reduced  wave 

z z 


equation: 


. k.2»  . 0. 

6x^  6y  ^ 


We  can  always  choose  the  amplitude  of  the  longitudinal  components  so 
that  one  of  the  transverse  components  is  zero.  The  choice  (H^  = 0)  cor- 
responds to  modes  with  polarization  predominantly  in  the  x-direction.  These 

modes  are  cabled  E modes.  The  choice  (H  = 0)  leads  to  modes  with  polari- 
pq  " y ^ 

zation  predominantly  in  the  y-direction.  These  modes  are  called  E^  , 

r J i pq  modes. 

In  this  paper,  we  consider  only  E modes,  but  the  sfune  analysis  can 
obviously  be  applied  for  E^^  modes. 


12 


! 


2.1.1  Waveguide  in  Region  (z  - 0) . 

Inside  the  core  we  can  write  the  electric  and  magnetic  fields  with  real 
propagation  constants  in  the  x-  and  /-directions.  For  the  two  logitudinal 


components  E and  H we  write: 
z z 


(13) 


(14) 


where  K and  K are  the  propagation  constants  in  x-  and  /-directions 

X / 


respective!/. 


This-  set  of  fields  corresponds  to  = 0.  For  the  reduced  wave  equa-: 


tion  (12)  to  be  satisfied,  the  propagation  constants  must  obe/  the  relation: 

(15) 

Substituting  in  (7)  to  (10),  we  find  the  transverse  field  components: 


„2_^2.2  .2  „2^^2 
K,  - n,  k - p_  - i\  + K 
11  Ox/ 


(16) 


(17) 


H.ho 


(18) 


H,= 


(19) 


Considering  relation  (3),  we  find  that  or  « 6^. 


If  we  compare  the 


electric  field  components  we  find: 


E % -§-  E 'V  — E . 

^ K,  ^ K^Ky  y 


(20) 


.1 

JL 


The  most  significant  field  component  is  E^,  which  means  that  the  wave  is 


polarized  predominantly  in  the  x-direction;  E^,  is  the  least  significant 


component  and  it  will  be  ignored. 

Outside  the  core  the  fields  must  vanish  at  a large  distance  from  the 
waveguide.  The  fields  must  be  represented  by  an  exponentially  decaying 
function,  which  means  that  the  propagation  constant  must  be  purely  imaginary 
We  assume  that  the  x-dependence  in  regions  3 and  5 is  the  same  as  in 
region  1,  and  that  the  y-dependence  in  regions  2 and  4 is  the  same  as  in 
region  1. 

With  the  requirement  that  the  tengential  electric  field  components  are 
continuous  across  the  core  boundaries,  the  field  components  outside  the 
core  can  be  written  as  follows: 

Region  2: 


(21) 


(22) 


"+fle  "]e  [e 

•1  p» 


(23) 


H.  5 0 


(24) 


with 


(25) 


(26) 


Region  4: 


r=^K(^K^)[e  ’+fle  "Je  ne  ] 

^kPa 

E,io  , HxSO 
^ ,Ur-t45- 


(27) 


(28) 


(29) 


with 


(30) 


In  regions  3 and  5,  we  adjust  the  amplitudes  so  that  the  strong  com- 
ponent E is  continuous  at  the  Core  boundaries.  E is  a second  order  quantity; 

X z 

thus  we  obtain: 


Region  5; 


I) 


(31) 


(32) 


(33) 
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Hy= 


"-^fle  ' ][e  -fle/ Je' 


(34) 


with 


(35) 


Y 


Region  5: 


E,  = -fl 


Cj  iO 


Hx  = o 


(36) 


(37) 


(38) 


(39) 


with 


Kj  = yi\i-P’=X^-'^s 


(40) 


We  note  here  that  the  E^  component  is  continuous  provided  that  we  neglect 
2 2 2 

K compared  to  n.  k in  the  numerator  of  equations  (33)  and  (3C). 

^ J 

We  know  that  the  tengential  field  components  must  be  matched  at  each 
boundary  separating  the  different  regions. 
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S SL 

In  the  x-direction,  at  the  boundaries  x = and  x = the  tangential 

field  components  are  E , H , E and  H . E is  already  matched.  E is 

z z y y z y 

negligible  and  by  matching  H , H automatically  will  be  matched. 

z y 

In  the  y-direction,  at  the  boundaries  y = -b  - and  y = -j,  the  tan- 
gential field  components  are  E , H , E and  H . E is  already  continuous, 

Z Z X XX 

is  zero,  E^  is  a second  order  quantity  and  we  will  neglect  it  here. 
is  the  only  remaining  tangential  component  to  be  matched. 

This  leads  to  the  following: 

^ a. 

at  X = -:r: 


(41) 


at  X = 


1^ 


. 'll?  - -lAfxf 


(42) 


a a 

X K "2  1 K 2 

Considering  e ^ and  e ^ as  unknowns,  eqs.  (41)  and  (42)  can  be  rewritten 
as  follows: 


/fx  <2 


(43) 


(44) 
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To  have  a solution  different  from  zero,  the  determinant  of  the  coefficients 
must  be  zero.  This  leads  to: 


-iT.  “''j 


(4S) 


Substituting  (45)  in  (41)  or  (42)  and  solving  for  tan  we  find; 


I y/  'n} ^*/3 


(46) 


This  is  the  eigenvalue  equation  of  the  TM  modes  of  the  infinite  slab 
waveguide.  In  the  particular  case  when  n.  = n_  we  have  y.  = y-,  and  sub- 

4f  fc  4f  ^ 

stituting  in  (45)  we  find: 

A"^  = 1 or  = ±1.  (47) 

For  A'  = -1,  we  substitute  in  (41)  and  we  find  the  even  modes  equation: 

a "l^  ^2 

tan  j = 2 “iT' 

n2  X 

For  A'  = +1,  we  have  the  odd  modes  equation: 

n ^ K 

<x  2 = (49) 

"l  2 

In  matching  in  the  y-direction  we  get: 


at  y = - (b  + y) : 


(50) 


at  y = 


-fl  . - a,[e“^S  ] 


(51) 


WSSSBBKEir^ 


From  (50)  we  find: 


H 


(C«l.)  _ ^ ^»^.ZC  a«w- 


From  (51)  we  find: 


c|^,C  » ^ ^.-3CWW 

Assuming  K"  ^ 0,  we  find: 


)(s 


(SI) 

(53) 
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(52) 


(S3) 


(54) 


or 


y 

2t|(^b  =2tancW^+2<:«^^^^+2L^ 


(55) 


which  is  the  eigenvalue  equation  of  the  TE  modes  of  the  infinite  slab  wave- 
guide. We  can  write  it  as: 

+ Yc) 

(56) 


For  the  particular  case  when  n^  = n^  we  have  Y3  = Yg  3^nd  (56)  becomes: 

2K  Yc 

tan  K b = — 

V 5 

or 


(57) 


tan  K 
y 

b y5 
2 " K“' 

(58) 

y 

K 

tan  K 

y 

1 

II 

n|<N 

(59) 

J 
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(58)  and  (59)  are  the  eigenvalue  equations  for  the  even  and  odd  TE  modes  of 
the  infinite  slab  waveguide,  respectively. 

From  now  on  in  our  analysis  we  will  assume  that  n2  = = n^  = n^, 

which  is  the  case  of  a rectangular  waveguide  immersed  in  a homogeneous  and 
isotropic  substrate. 

The  constant  k"  can  be  determined  as  follows. 

Let  us  first  consider  the  even  modes  in  the  y-direction.  We  can  write  (58)  as; 


Yc 


^-2i  arctan  _ -i  Ky^. 

' y ~ ® 


(60) 


substituting  in  (53)  we  find: 


For  the  odd  modes  we  have: 


^ i^y(c+b) 


2i  arctan  -i  K„b 

e Yr  = e y ; 


(61) 


(62) 


also,  we  can  write  (53)  as: 


i KyrC  2i  arctan 

e y = -A  e 


(63) 


substituting  (62)  in  (63)  we  find: 

(64) 

Finally,  we  see  that  the  constants  k'  and  k"  can  be  determined  from  (47), 
(61)  and  (64)  depending  on  the  kind  of  modes.  The  only  remaining  unknown 
constant  in  the  expressions  of  field  components  is  A,  which  will  be  deter- 
mined from  the  amount  of  power  traveling  in  the  waveguide. 

The  eigenvalues  for  the  propagation  constant  determined  by 
(46)  and  the  eigenvalues  for  determined  by  (57)  remain  the  same 
along  the  z-axis  as  long  as  there  is  no  change  in  the  physical  parameters 
of  the  waveguide. 
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or  graphically.  Using  the  second  way,  the  solution  for  the  eigenvalue 

equation  (57)  is  presented  in  Fig.  4.  As  a variable  we  use  x = Kb,  and  the 

' 2K  Yr 
y ' 5 

solution  is  determined  by  the  intersection  of  the  branch  tan  x with  — — =• ; 

K -Yc 
y ' 5 

Yg(x)  can  be  obtained  by  subtracting  eq.  (25)  from  eq.  (15): 

<•2  2.,  2 i,  2 2 

(nj  - n^  )k  = Ky.  + Yg  , (65) 

or: 


and: 

1^3  ' 


(66) 


(67) 


This  solution  determines  the  propagation  constants  in  regions  (z  5 0) 

2 2 2 2 

and  (z  - L) . Some  of  the  curves  in  Fig.  4 were  drawn  for  (n^^  - U2  )k  b = 100. 

Under  this  condition  the  modes  are  cut  off  for: 

=/l00^’  ^ 0 

or: 

X ^ 10.  (69) 

We  see  that  there  are  eight  intersections,  and  this  guide  can  carry 
eight  modes. 

2 2 2 2 

For  (n^  - n^  )k  b =10,  the  curves  drawn  in  Fig.  4 show  that  only  four 

2 2 2 2 

modes  can  be  carried  by  the  guide.  By  decreasing  (n^  - n^  )k  b , we  decrease 

the  number  of  possible  modes  in  the  guide. 


; 

1 

I 
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When  we  reach  the  region  (z  - 0)  (see  Fig,  3),  we  have  both  guide  I and 
guide  II  present.  '!"he  boundary  conditions  in  the  x-direction  are  the  same 
as  in  region  (z  - 0);  as  a consequence,  the  modes  in  the  x-direction  and 
the  eigenvalues  for  will  remain  the  same  as  in  region  (z  - 0).  In  the 
y-direction  we  will  have  an  exchange  of  power  between  the  two  waveguides 
through  the  separation  region,  by  means  of  evanescent  waves.  The  boundary 
conditions  which  we  must  impose  on  the  fields  in  this  region  are  different 
from  the  conditions  we  imposed  in  region  (z  - o) ; as  a consequence,  the 
eigenvalue  equations  (58)  and  (59)  are  not  valid  in  this  region,  and  the 
propagation  constant  will  be  determined  by  new  eigenvalue  equations.  The 
new  Ky  will  be  a function  of  the  separation  distance  c.  We  note  here  that 
the  propagation  constants  in  region  (z  ^ L)  are  the  same  as  in  region  (z  - 0) . 

The  above  analysis  shows  that  light  will  propagate  in  region  (0  ^ z ^ L) 
with  a propagation  constant  different  from  that  of  regions  (z  - 0)  and  (z  - L) . 
This  means  that  we  will  have  a reflection  in  guide  I at  :he  boundary  (z  = 0), 
and  a reflection  in  guide  II  at  the  boundary  (z  = L).  Another  effect  can 
take  place  if  the  guides  were  multimode  carriers.  This  effect  is  the  exci- 
tation of  undesired  modes  and  coupling  between  different  modes. 

For  maximum  power  transfer  from  guide  I to  guide  II,  the  two  reflections 
at  (z  = 0)  and  at  (z  = L)  must  be  minimized,  which  means  that  we  have  to  look 
for  optimal  conditions  so  that  the  difference  between  the  propagation  con- 
stants in  the  different  regions  along  z is  minimum.  Also,  the  excitation 
of  undesired  modes  can  be  eliminated  if  the  propagation  constants  of  the 
desired  modes  are  the  same  in  all  three  regions. 

In  the  light  of  the  above  discussion,  we  need  to  find  the  propagating 
modes  in  the  three  different  regions.  This  is  done  by  solving  the  corre- 
sponding eigenvalue  equations.  The  solution  can  be  obtained  either  numerically 


“■m 


In  the  region  (0  - z - L),  the  solution  will  be  presented  in  Chapter  HI 


2.1.2  Waveguide  in  Region  (z  - L) . 

For  guide  II  in  region  (z  - L)  and  for  both  guides  in  region  (0  - z - L), 
the  boundary  conditions  in  the  x-direction  are  the  same  as  for  guide  I in 
region  (0  - z).  As  a consequence,  the  x-dependence  of  the  field  components, 
the  eigenvalue  equations  for  modes  and  the  propagation  constants  in  the  x- 
direction  are  the  same  in  all  three  regions.  From  now  on  we  drop  the  x- 
dependence  of  the  field  components  in  our  equations. 

However,  in  the  y-direction  the  boundaries  are  not  the  same  in  all  three 
regions,  and  therefore  the  y-dependence,  the  eigenvalue  equations  for  the 
modes  and  the  propagation  constants  also  vary  from  region  to  region. 

For  the  case  of  guide  II  in  region  (z  ^ L) , we  can  write  the  field 
components  as: 

Region  1: 


, H.30 


• 7 •- 


1 
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Region  3: 


_2  r » tK,(b*i),  -iTjCa-b-f) 

e,=8t[c  +8fi’  Me  (74) 

* m| 


* ^ Rjf 

£,iro  , H.=o 


(76) 


(77) 


Region  5: 


E*  ='6|f^‘  Ce J e 

E,S(0  , H,  =0 


H is  the  only  component  which  has  to  be  matched  at  y » -s-  and  at 
z 2 

c c 

y*b  + -2.  Aty  = b + -^we  find: 


(78) 


(79) 


(80) 


A 


i 


t 

t 

I 

i 

t 

j 

j 


(81) 


or: 


At  y = -j  '*'®  find: 


* kWUw  ^ 

e ’ =8e  (i, 


Assuming  Q"  ^ 0 and  dividing  eq.  (85)  by  eq.  (83): 

Ztl^kb  7^  tliOAcloM.^ 

e =e  «i  =<2  K) 


2 (6,  2 ■ a-S 

Eq.  (87)  is  the  same  as  eq.  (57),  which  means  that  in  the  y-direction  the 
modes  of  guide  II  in  region  (z  - L)  are  the  same  as  the  modes  of  guide  I 
in  region  (0  - z). 

The  constant  B"  can  be  obtained  from  (85);  it  is  given  by: 


B"  = e 


•iKy  (b+c) 


• j::-?  ■ 
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for  even  modes,  and  by; 


for  odd  modes. 


2.2  Power  in  Rectangular  Waveguide 

If  we  are  far  from  cut-off,  the  fields  outside  the  core  are  very  small 
compared  to  the  field  inside  the  core,  and  almost  all  the  power  is  contained 
inside  the  core.  As  an  approximation,  we  will  assume  that  all  the  power  is 
contained  inside  the  core. 

The  power  can  be  obtained  by  integrating  the  z-component  of  the  power 
flow  vector  (Poynting's  vector) 

= I Re(E  X H)  • e^  (91) 

over  the  transverse  section  of  the  waveguide  core.  Eq.  (91)  implies  that 
is  a time-averaged  qua'ntity. 

In  our  analysis,  we  considered  no  reflected  wave  in  the  z-direction, 
and  Bq  > 0.  If  Pj  is  the  power  in  guide  I,  we  thus  obtain: 


Pj  = if  ] , 

-2"  -b-i 


I 2 I 

but  we  know  that  H = 0 and  E = 0,  and  with  the  help  of  (16)  and  (19)  we  find: 
X y 

I C • ® Pa  Ktf  J 


-U-i- 
1 b 2 


I 
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u. 

fJ=-i|fliyj«,^[-«rk-K.jJ  [-i-lfl'I+fle  +rt'e  JctX 


Ci 

■ 

2 


(94) 


We  consider  relations  (47),  (61)  and  (64),  the  constant  A'  = -1  or  +1  cor- 

j A"  iKv(b+c)  -ifCv(b+c) 

responding  to  even  or  odd  modes  respectively,  and  A =e/'  or-e  ^ 
corresponding  to  even  or  odd  modes  respectively,  and  we  substitute  A'  and  A" 
in  (94) . We  find  that 

Ot 


c 


o. 

f 


-b-l 


After  performing  the  integrations  we  find: 


K»j 


or: 


(95) 


(96) 


(97) 


The  upper  sign  corresponds  to  even  modes  and  the  lower  sign  corresponds 
to  odd  modes.  Substituting  for  tan  — =—  and  tan  their  respective  values 


from  (48)  or  (49)  and  (58)  or  (59),  the  right-hand  side  of  (71)  is  unchanged, 
and  we  finally  have  that: 


Pi 


(98) 


By  solving  Maxwell's  equations  and  imposing  the  boundary  conditions, 
we  were  able  to  describe  the  electric  and  magnetic  fields  in  all  regions, 
and  to  determine  the  values  of  all  constants  and  parameters  except  for  the 
constant  A.  The  magnitude  of  this  constant  can  be  determined  from  eq.  (98) 
when  we  know  the  amount  of  power  traveling  in  the  waveguide.  We  see  that 
if  Lhe  propagation  constants  remain  constant  along  the  waveguide,  the  power 
is  proportional  to  the  square  of  the  absolute  value  of  A. 

In  the  region  (z  - L),  the  constant  B is  related  to  the  power  inside 
guide  II  by: 


J 


(99) 


2.3  Reflection  Coefficients. 

An  important  parameter  which  will  be  used  frequently  in  the  following 
is  the  reflection  coefficient  at  the  boundary  (y  = for  guide  I and  at 
the  boundary  (y  = -j)  guide  II. 

To  compute  this  parameter  in  the  case  of  guide  I in  region  (z  - 0),  we 
see  that  the  y-dependence  of  the  electric  field  (see  relations  (13)  and  (16)) 
is  given  by: 


I 


(100) 


Considering  the  first  term  as  incident  field  and  the  second  term  as 
reflected  field  E , we  see  that  the  reflection  coefficient  at  the  boundary 


(y  * --j)  is: 


t;  lu--^ 


(101) 


where  the  upper  sign  corresponds  to  even  modes  and  the  lower  sign  to  odd  modes. 

In  guide  II,  the  incident  and  reflected  electric  fields  at  the  boundary 
(y  = and  e respectively. 

The  reflection  coefficient  at  this  boundary  is  given  by:  ' 


Ri  = rl  c =-rj 


C 3=2 


n te 


(102) 


where  the  upper  isgn  corresponds  to  even  modes  and  the  lower  sign  to  odd  modes, 
We  see  that  Rj  = Rjj  = R;  this  result  is  due  to  the  fact  that  guides  I and  II 
are  symmetric.  Also,  we  observe  that: 


7 i |2 

\R\^-.\te  ’ 1 


(103) 


which  indicates  a total  internal  reflection. 

In  the  region  (0  - z - L)  the  situation  is  totally  different.  Here  both 
waveguides  are  present  and  there  is  power  coupling  between  them;  the  reflec- 
tion coefficients  at  the  boundaries  (y  = ")  and  (y  = j)  are  smaller  than  R 
in  absolute  value.  To  improve  the  boundary  conditions,  we  have  to  consider 


J 


CHAPTER  III 


TWO  PARALLEL  RECTANGULAR  WAVEGUIDES 


In  this  chapter  we  will  derive  the  eigenvalue  equation  for  modes  in 
the  region  (0  - z - L) . We  will  then  present  the  method  of  Miller  and 
Marcatili  for  the  solution  of  coupling  between  the  two  waveguides  in  this 
region,  and  discuss  the  validity  of  their  approximations. 

The  reflection  and  transmission  coefficients  at  the  boundaries  (y  = -y) 
and  (y  = -j)  ^re  important  parameters  to  find  the  eigenvalue  equations  for 
modes  and  to  find  the  coupling  between  the  two  guides. 

In  the  case  of  the  single  rectangular  waveguide  studied  in  Chapter  II, 
let  us  take  a close  look  at  the  eigenvalue  equation  (55).  We  see  that  the 
condition  for  guided  modes  can  be  stated  as  follows: 

2iK  b = 2iirn,+  the  phase  of  the  reflection  coefficient  at 

' c 

the  boundary  (y  = -j) 

+ the  phase  of  the  reflection  coefficient  at 
the  boundary  (y  = j) . 

To  find  the  eigenvalue  equation  for  modes  for  guides  1 and  II  in  the 


region  (0  - z - L) , we  will  use  the  same  method  already  used  in  Chapter  II. 

For  guide  I the  reflection  coefficient  at  the  boundary  (y  = -j  - b)  is  the 
same  as  in  region  (z  - 0),  but  because  of  the  presence  of  guide  II  the  reflection 
coefficient  at  the  boundary  (y  = is  different  from  its  value  in  the 
region  (z  ^ 0). 

3.1  Reflection  and  Transmission  Coefficients  in  the  Region  fO  ^ z ^ L)  at 


the  Inner  Boundaries 


■^fLand.iY  =.|I. 


The  two  waveguides  are  identical,  so  the  reflection  coefficient  in 
guide  I at  the  boundary  (y  = is  the  same  as  the  reflection  coefficient 
in  guide  II  at  the  boundary  (y  = -j) . 


BnBSS 
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at 


y - y: 


£2  e' 


Therefore,  by  dividing  eq.  (113)  by  (112): 


/ c- 

^5 


(112) 

(113) 


(114) 


From  (114)  we  find  that: 


I..'  I^'a 


l+t  ia 

ys 


Kb 


(115) 


X54*»'I^3 


and,  therefore,  from  (110)  and  (111): 

Ec-Ey . !dL^  ^ - c li.  .— 
E;*E,  ■ Kj  *^1  /*r 


(116) 


Cr 

The  quantity  ^ is  the  reflection  coefficient;  from  (116)  we  find: 


M 

l + R' 


o'  J / L 


.tc  ys+cKh-cys”*’^^ 
r-*-  — - ^ — t: — rr  . -i¥tC 


*■  J 

..  Xs  IjSWiJsC+tlCjCkJfsC 


(117) 


<5  Vst*'-'’  KiJSU)ist 


y5cm5^'*’*''K5  (118) 


(119) 


- - 
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or: 


I'Kj  K^CXx  fc5  C 


R c 


K}-is  ^ ihl^  . |i  V 

fiMt 


(120) 


(121) 


- ^ V /o"  he)  . — 


We  put: 


II  - SU^sC- 

■ VCKs^Ks^) 


and  finally  we  find: 


1(4  Is 


/i+a? 


The  transmission  coefficient  is:  t = — . 

Ei 

From  (112)  and  (113)  we  find: 


c z&e***^ 

t =‘,  . (Cs, 


(123) 


(124) 


(125) 


from  (110)  and  (111)  we  have: 


From  (125)  and  (126)  and  using  -g-  = R'  we  find: 


l*i 

^5 

Substituting  R'  from  (120)  in  (127)  we  find: 


i = e 


f <.Us +K»‘l  SU!!sC  (.Kj-i  Its) 


t = e ***|  [lXjlCiC(,t5C.  - K^(KjV  lls'j  Skl(5« +»?1  (Ks'+Vsh  SKI!  jc] 
t[l(s(Kj^^SKl!sCf2lfs<<ii‘ai(5C+J5(«;5VI(s';sKI!^  j 

|^21!5l^CUT!jC+t  ^^54. 1 1^3)*^ 

t : £^‘'[2K^!s  (CKIst+Sk)(sCj  4.uK5Kj*(cU)!jC+SU5cjl  X 
^(l!5+tl!^}|lVf([^CliVsC'  + <'CK^^^|)5!i 


1 
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(126) 


(127) 


(128) 


(129) 


(130) 


Now  substitute  chy-c  + shYpC  = e in  (130);  we  find 


J^%l<l3CUK5C4.tCl(s-»?')SW>!isC 


Comparing  (131)  with  (120),  we  find 


so  that  energy  is  conserved 


We  now  consider  two  limiting  cases 


a)  For  infinitely  large  separation  distance  we  have 


u = oo  and  cotanh  Yc^  = 1 


For  zero  separation  distance  (c  = 0) , the  tvfo  guides  touch  each 
other  and  form  one  guide  with  a width  equal  to  twice  the  width  of 
the  original  guide.  In  this  case  we  have: 
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substituting  in  (124)  and  (125),  we  find; 


R'  = 0 


and: 


/ i‘cVl£taM>C  , 


= -1 


(136) 


An  important  remark:  We  see  that  at  any  given  distance  c we  have 

a phase  difference  of  it/ 2 between  reflection  and  transmission 
coefficients. 


3.2  Guided  Modes  of  the  System  of  Two  Guides  in  Region  (0  ^ z ^ L) . 
The  condition  for  guided  modes  in  guide  I is  the  following: 


2 i K b 

y 


the  phase  of  the  reflection  coefficient  at  the.  boundary 
(y  = -b  - — ) + the  phase  of  the  reflection  coefficient 
at  the  boundary  (y  = -j)  + 2imr 

phase  of  R + phase  of  R'  + 2imr 


V-*5  '' 


or: 


The  solution  of  this  eigenvalue  equation  determines  the  propagation 
constants  of  the  propagating  modes  in  the  region  (0  ^ z ^ L) . In  Fig.  S, 
we  plotted  tan  2Kyb  and  F(Kyb)  in  function  of  the  variable  x = K^b  and  for 


(137) 


(138) 


0.1  T) 
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different  values  of  the  distance  c.  Each  intersection  of  the  two  above  func- 
tions determines  a mode  and  a propagation  constant  K^. 

2 2 2 2 

We  see  that  for  (n^  - n2  )k  b = 100  we  have  eight  modes,  and  for 


Cn, 


2 2 2 2 

- n-  )k  d'^  = 10  we  have  four  modes.  The  modes  are  cutoff  when  K - 


k/n^^  - By  comparing  Fig.  4 with  Fig  5,  we  see  that  the  propagation 

constant  varies  very  slowly  with  the  separation  distance  c.  This  means 
that,  except  when  c is  quite  small,  an  incident  mode  propagates  approximately 
unchanged  into  the  coupling  region.  In  other  words,  the  coupling  coefficients 
of  the  incident  mode  to  the  modes  that  can  exist  in  the  coupling  region  are 
extremely  small  except  for  the  mode  whose  is  almost  identical  to  the  of 


[ the  incident  mode  (in  this  case,  the  coupling  coefficient  is  very  close  to 

i-  2 2 2 2 

i unity).  For  example,  if  (n^^  - n^  )k  b = 100  or  10,  then  the  above  statements 

I are  true  for  all  c - b/4  or  c - b,  respectively. 

In  the  region  (0  - z - L) , we  have  exchange  of  power  between  the  two 

guides,  which  means  that  the  power  in  each  guide  is  a function  of  z.  With 

the  assumption  that  the  x-  and  y-dependence  along  the  z-axis  does  not  change, 

the  power  in  each  guide  is  proportional  to  the  square  of  the  absolute  value  of 

2 

the  field  amplitude.  The  power  will  be  proportional  to  |a|  and  will 

2 

be  proportional  to  )b)  . Our  problem  is  to  find  the  z-dependence  of  the 
I amplitudes  A and  B. 

I Dropping  the  x-  and  y-dependence,  the  electric  field  can  be  written  as: 

^ -iSz 

I Ej  = A(z)  e , in  guide  I,  (139) 

E2  = B(z)  e in  guide  II,  (140) 

where  B is  the  propagation  constant  in  the  z-direction  in  region  (0  - z - L) . 


I 


I 


9 

« 

J 

I 
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3.3  Miller-Marcatili  Solution. 

S.  E.  Miller  [2]  studied  the  coupling  between  transmission  lines.  His 
solution  was  used  by  Marcatili  [3]  and  applied  to  the  coupling  between 
rectangular  waveguides. 


1 


Miller  set  up  two  differential  equations  for  the  field  amplitudes  in 
guide  1 and  guide  II: 


6^1  6^2 

where  and  are  the  rates  of  change  in  the  electric  field  in  guide  I 
and  guide  II. 

If  the  two  waveguides  are  identical  we  have  that  \~i  ~ \~2  ~ I “ iS, 
the  propagation  constant  in  the  z-direction  in  the  coupling  region;  = k^2 
represents  the  reaction  of  the  coupling  mechanism  on  the  guide  itself,  and 
k22  = = k represents  the  transfer  effect*  of  the  coupling  mechanism. 

Miller  assumes  a very  large  separation  distance  between  the  lines;  this 
means,  in  our  case,  that  the  cliange  in  the  propagatio-  constants  between  the 
regions  (0  - z)  and  (0  - z - L)  is  negligible  and  the  coupling  between  the 
two  guides  is  very  weak  so  that  no  reflected  wave  xs  excited.  Also,  Miller 
assumes  that  k = kj^^  “ ^21  ^ purely  imaginary  quantity,  which  is  a 

necessary  condition  to  have  a solution  to  equations  (1')  and  [2')  in  the  form: 


= cos  Cz  e 


-(iC+Y)z 


-(iC+Y)z 


E2  = i sin  Cz  e ^ . (4') 

Here  E^  and  E^  are  normalized.  C = -ik  is  the  coupling  coefficient,  and 
0Q  “ -iy  is  the  propagation  constant  in  the  z-direction  in  the  region  (0  - z) 
Marcatili  [3]  considers  C as  the  change  in  the  propagation  constant  in 
the  z-direction  between  the  two  regions  CO  - z)  and  (0  - z - LJ , and  he 
calculates  C by  solving  Maxwell's  equations  in  the  region  (0  - z - L)  for 
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the  system  of  two  parallel  waveguides.  By  imposing  the  boundary  conditions, 
he  finds  the  propagation  constant  6,  then: 


/•  rt  « . I 2.1^3  05  C 

C c S-Po  — rr — 3 • 

' P«b 


The  value  of  the  distance  L necessary  for  a complete  transfer  of  power  from 
guide  I to  guide  II  corresponds  to: 

lE2l  = l=sinI  (6') 


I . CL  and  L . i - ’8b  «/  > y/) 

"VS 

Marcatili  points  out  that  his  solution  is  valid  only  for  a very  large 
separation  distance  c,  for  which  the  distance  L given  in  {!')  equals  several 
hundred  wavelengths. 

We  want  to  investigate  the  validity  of  Miller's  solution.  He  implies 
that  k is  a purely  imaginary  quantity;  if  we  investigate  this  quantity 
closely,  we  see  that  it  is  given  by  the  relation: 

where  t is  the  transmission  coefficient  and  <J>  is  the  phase  change  when  the 
optical  ray  travels  the  path  between  two  successive  incidences.  4>  cannot 
be  known  unless  we  know  the  z-dependence  of  the  field  components. 

In  our  method,  we  use  real  quantities  in  setting  up  differential  equa- 
tions, and  these  are  known  without  imposing  a type  of  solution  in  advance. 

We  show  that  Miller's  assumption  that  k is  purely  imaginary  is  not  valid. 

(Xir  solution  will  remain  valid  for  a small  separation  distance,  and  there- 
fore strong  coupling,  between  the  guides. 


CHAPTER  IV 


COUPLING  BETWEEN  RECTANGULAR  WAVEGUIDES 


4 . 1 Coupling  Mechanism. 

In  Fig.  3 we  represent  the  system  of  two  parallel  identical  rectangular 
waveguides  in  the  y-z  plane.  We  assume  the  two  guides  to  be  extended  from 
(z  = -®)  to  (z  = +“),  with  only  one  excited  mode  traveling  in  the  positive 
z-direction  in  both  guides. 

In  guide  I we  will  follow  the  path  of  an  optical  ray  with  associated 
electric  field  E^.  This  ray  is  partially  reflected  at  the  interface  (y  = 
at  the  point  z = z^  and  partially  transmitted  to  guide  II.  Also  at  the  same 
point  z,  in  guide  II  we  have  an  incident  electric  field  E2  at  the  interface 
(y  = + -2)-  We  can  write  the  two  fields  and  E2  at  the  point  Zy  in  terms 
of  Ej^  and  E2  at  the  same  point  as  follows: 


Ej"  = R"Ej  + tE2 


E2'  = R'E2  ^ tE^ 


(141) 


(142) 


R'  and  t are  the  reflection  and  transmission  coefficients  at  the  boundaries 
c c 

(y  = -rp  and  (y  = j) ; they  are  the  same  for  both  guides . 

The  fields  and  E2'  are  totally  reflected  at  the  boundaries  (y  = -(b  "j)) 

and  (y  ® b * p»  then  they  meet  again  at  the  coordinate  Z2;  at  this  point 
the  two  fields  can  be  written  as: 


(143) 


(144) 


J 
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where  4>  is  a constant  phase.  and  are  the  incident  electric  fields  in 


guide  I and  guide  II  respectively  at  the  point  z^,  while  Ej^' 


z=z. 


and 


'2  'z=z. 


are  the  fields  obtained  after  traveling  one  bounce  along  the  z-axis. 


We  have  that 

Az  = Z2  - z^  = (145) 

Firstly,  we  impose  energy  conservation.  We  know  that  the  power  in 

2 

guide  I is  proportional  to  |Ej^|  and  the  power  in  guide  II  is  proportional 
2 

to  IE2I  . With  Ej  and  E^  normalized,  the  total  power  will  be: 


’'0  "'l  * ‘’ll  ’ •. 

We  want  to  verify  that  [E^  * 1^2  " I'^l  * 1^2  ^^1^ 


We  can  always  write: 


E2  = a e^®E 


(146) 

i2 


(147) 


with  a(z)  and  9 (z)  functions  of  z. 

Substituting  (147)  in  (141)  and  (142),  we  find: 


Ej  (Z2)  = (R"E^  + tae^®  E^  e^® 


i9 


i0. 


, e (R  + t a e ) . 


and 


E2  = (R'a  e^®  Ej  + t Ej)  e^®  = E^  e^®  (R"a  e^®  + t) . 


lE^  (Z2)l^  + IE2  (Z2)|^  = |Ej^  (|R'  + tae^®!^  + [R'a  e^®  + t[^)  (150) 


(148) 

(149) 


[ja+i-ae  ®|V  loae'®-c|^  j 

r'ilf  +a?u?  cje  + ] 

i+u^  L 
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Finally, 


I " 

Po(^)=  P.W 


Thus,  we  see  that  the  total  power  is  constant  along  z. 


4 . 2 Differential  Equations  for  Power  Transfer. 

In  calculating  the  power  in  each  guide  as  a function  of  z,  we  deal  with 
the  absolute  values  of  the  field  components  and  can  thus  eliminate  some  of 
the  phases  involved,  particularly  the  phase  <p  which  is  as  yet  unknown. 

We  assume  that  the  bounce  length  Az  = Z2  - Zj  is  very  small  compared  to 
the  distance  L,  and  that  the  separation  distance  c is  large  enough  so  that 
we  can  consider  the  power  exchange  per  bounce  as  a very  small  quantity.  Also, 
we  assume  that  the  propagation  constants  in  the  three  regions  (0  - z) , 

(0  - z - L)  and  (z  - L)  are  the  same,  so  that  there  is  no  reflection  at  the 
boundary  (z  = 0)  in  guide  II,  and  we  neglect  the  reflected  wave  in  guide  I 
at  the  boundary  (z  = L) . This  last  assumption  means  that  either  the  power 
Pj  is  negligible  at  Cz  = L),  or  guide  I is  tapered  about  (z  = L) , or  both. 

The  field  amplitudes  in  guide  I and  guide  II  are  called  A and  B respec- 
tively, as  in  Chapter  II.  After  dropping  the  x-  and  y-dependences,  which 
are  the  same  as  in  region  (z  - 0)  for  guide  I and  as  in  region  (z  - L)  for 
guide  II,  the  fields  in  region  (0  ^ z ^ L)  are  given  by: 


Ej  = A(z)  e 


(151) 


= B(z)  e 


(152) 


A and  B are  functions  of  z because  of  the  coupling  between  the  two  guides. 
Here  we  consider  the  normalized  value  of  A and  B and  we  have  the  following 
relations: 


I 
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Pj(z)  = lA(z)|^ 
P,-(z)  = |B(z)l^ 


The  change  of  Pj  per  unit  length  along  z is  given  by: 


45 . !?(?.)-?(»■). 


dz  Za  - Zi 


(153) 


We  can  write: 


B(z)  = a(z)  A(z) 


substitute  (155),  (124)  and  (132)  in  (154),  and  find: 

jz  *pb'  '''*  j 


. j4 P/.,+  a*+2U.a*^«)=  Dr-R+CL*l5*2u.a>«-»?] 

■ 2pb(H-li*)  ‘ ^ '• 

(1 


r 


p. 


II  ' fB 
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where  the  constant  D = 


-j-  . For  Pjj  we  have: 


2eb(l+u  ) 

^ (la(Ji)  l - 1 atz,  » A [^1 6 1?'+  -t  fl  1^- 1 B I*-  J 

1/  l\ 

= D Ifl  > I J 

. P,fl|ZraWe«Je+  u.VA,le  + 2u.o-  »«9+|-a?  U+a*)j 
«(z  ' 


(161) 


(162) 


We  see  that  ■■  + = 0,  so  that  energy  is  conserved.  To  solve  the  two 
differential  equations  (160)  and  (162),  we  must  find  the  function  sin  0(z), 
where  9(z)  is  the  phase  difference  between  B and  A.  IVe  write: 


reczo-ecz,)]  . 
dz  zfih  i j ’ 


(163) 


0(Zj^)  is  the  phase  difference  between  B(Zj^)  and  A(Zj),  and  0(Z2)  is  the  phase 


difference  between  B(z2)  and  A(z2) • 


g (zo  = phase  o|  6 (^i ) - ^ 

spls«erR'6(z,)+-trtC*iJ]e  - ph<we(R'if!«)+l6t*'>]e 

= pW (o-e  - ^ ) e ^ ^ _ 


/rt  (9  c ouue.f g**. ( L!^i 

_ phftieO-  •' 


a«.»_/Sg!g1 


ecz,)  = a«t«.^ 


(164) 


J 


therefore: 


, ra/>v»v^~ 

■tovw  0CZi)  = Latflse 


J/  ^ 


u?-l  . gd-a-V  _!_ 

= S?:7  acavi;  «se 


(165) 


and 


dltctM^  ^iAw.9  /-2t<jA.v9-U- ). 

2j3bC+‘^")  ^ ^ 


olz  tz  Zi  - z 

Finally,  substituting  D = — 


(166) 


K Pjj  2 

^ and  ■= — = a , we  find: 

26b  (W)  *^1 


PCCgfl(2wg-^^y~^  ) 

6<Z 


(167) 


Eqs.  (160),  (162)  and  (167)  are  three  differential  equations,  which  are  to 
be  solved  simultaneously.  These  equations  represent  our  general  formulation 
of  the  coupling  between  two  waveguides. 

4 . 3 Particular  Cases. 

We  discuss  two  particular  cases. 

a)  Assume  that  at  an  arbitrary  point  z it  happens  that  the  fields  in 
the  two  waveguides  are  the  same  in  phase  and  magnitude;  then  A = B. 
At  the  next  incidence,  i.e.  after  one  bounce,  we  have  that 
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A(z  + Az)  = [A(z)R^  + A(z)t] 
B(z  + Az)  = [A(z)R'  + A(z)t] 


(169) 


The  two  fields  have  the  same  phase  and  same  magnitude,  and  the  waves 
propagate  in  the  two  waveguides  as  if  there  were  no  coupling 
between  them. 

b)  This  case  applies  to  the  geometry  of  Fig.  3.  We  assume  that  at  a 
given  point  z in  the  field  amplitude  B is  leading  A by  then  at 


the  point  z + Az  we  have: 


A(z  + Az)  = [R^A(z)  + tB(z)]  e , 
B(z  + Az)  = [R'B(z)  + tA(z)]  e^*^. 


(170) 

(171) 


e cr.+w> = [Ri.rtc*)e  4 e ' *-;R/K«e  Va) 


(172J 


(173) 


From  (172)  and  (173)  we  see  that  the  magnitude  of  A decreases,  and 

the  magnitude  of  B increases  but  the  phase  difference  remains  the 
37T 

same,  i.e.  ^ , which  means  that  the  phase  6(z)  remains  constant 


along  z,  and 


sin  0 = -1. 


(174) 


An  important  observation  is  that  when  the  transmitted  amplitude 
from  guide  I to  guide  II  and  the  reflected  amplitude  in  guide  II  are 
in  phase,  the  total  amplitude  is  the  sum  of  their  absolute  values. 
But  if  the  transmitted  amplitude  from  guide  II  to  guide  I and  the 
reflected  amplitude  in  guide  I have  a phase  difference  of  u,  then 
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the  resulting  total  amplitude  is  the  difference  between  their 
absolute  values. 

Substituting  (174)  in  (160)  and  (162),  we  find  the  differential 
equations  for  the  coupling  between  the  two  rectangular  waveguides 
in  the  form: 


<iPl 

dZ 

dlPp 

dz 


rD[Pj-PB+ZU.|/f^] 


(175) 


(176) 


These  two  differential  equations  are  valid  as  long  as  the  phase  of 

the  field  in  guide  II  is  leading  the  phase  of  the  field  in  guide  I 
3tt 

by  ~Y.  However,  when  the  magnitude  of  A is  so  small  that  the  trans- 
mitted field  from  guide  II  to  guide  I is  larger  in  magnitude  than 
the  reflected  field  in  guide  I,  then  the  phase  of  the  resultant 
field  in  guide  I will  lead  the  phase  of  the  field  in  guide  II.  In 
this  later  situation,  the  power  begins  to  flow  back  from  guide  II 
to  guide  I,  and  eq.  (174)  becomes: 

sin  6 = +1  (177) 

Substituting  (177)  in  (160)  and  (162),  we  get  the  differential 
equations 


dZ 

^-D[P„-Pi42U^J, 


(179) 
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which  differ  from  (175)  and  (176). 


The  change  in  the  phase  9 from  to  occurs  when: 


|r'A|  ^ |tB|  or  |R"a(^  ^ |tB(^  or  |R"1^Pj  ^ 


(180) 


Substituting  (146)  in  (180)  we  get: 

— iPi  <-^:(Po-Pl) 


(181) 


or: 


p y , and 

I y i+u.*’  i+u’’ 


P 

P,4 

1+u? 


(182) 


4.4  Solution  for  Power  Transfer. 


4.4.1  Approximate  Solution. 

In  these  two  differential  equations,  let  us  compare  Pj 
with  the  first  bounce  after  z = 0.  We  have: 


Pjj  to  2u/PjPjj, 


Pl=lR'fPo=4Po  , 


p « To 

In  — 


/+u> 


(183) 


and 


fi-Pn=—  > 

® u.V>  ^ ^ i+u?-  i+a*- 


(184) 


We  see  that  Pj  - Pjj  < h(.2\xy/P~^~^)  and,  as  the  wave  propagates  along  z,  Pj 


decreases,  Pjj  increases,  Pj  - Pjj  decreases  and  2u/PjPj j increases.  If  we 
neglect  Pj  - Pjj  compared  to  2u/PjP j j in  (175)  and  (176),  the  two  equations 


become: 


4 Pi 


d2 

4 Pi 


= -2uj>y^^ 


(185) 


dz 


(186) 


"1 


f 
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With  the  initial  condition  Pj  = and  P^^  = 0 at  z = 0,  the  solution  of  (185) 
and  (186)  is: 


Pj  = Pq  cos  uDz, 

2 

Pjj  = Pq  sin  uDz; 


(187) 

(188) 


this  corresponds  to: 


|A(z)|  = ^ cos  u Dz, 


or  A(z)  = A(0)  cos  u Dz, 


and,  remembering  that  B leads  A by  , we  have: 

B(z)  = -iA(0)  sin  u Dz, 

K K ^ 

y u y ' 5 

The  factor  uD  = x tt,  with  u = ' shYrCJ  for  a sufficient  large 

® Vs  ?jC 

separation  distance  we  have  u -►  ■»,  and  therefore  U.  d '■■■  ^ and 

216,  <5 


(189) 


(190) 


2pb  U.  (KjVjj')PW 


thus,  we  see  that  this  factor  is  the  same  coupling  coefficient  C of  eq.  (5'). 
Finally,  we  observe  that  the  approximate  solution  (189)  - (190)  is  the 
Miller-Marcatili  solution,  and  that  it  is  obtained  by  neglecting  the  term 
Pj  - Pjj  in  our  differential  equations.  Assuming  that  the  factor  k used  by 
Miller  is  a purely  imaginary  quantity,  and  observing  that  setting  1 - R'  = t 
is  the  same  as  neglecting  the  term  Pj  - Pjj,  we  conclude  that  Miller's 
a-proximation  is  not  accurate  except  for  a very  large  separation  distance 
between  the  two  guides,  in  which  case  Miller's  treatment  approaches  our  more 
rigorous  formulation. 


4.4.2  Exact  Solution. 


We  introduce  a new  function  Z(z)  given  by: 


then: 


wx  • 

<lz  ^*Z  Vp,  dz 

Using  (176),  (191)  and  (192),  we  find: 


or 


zZdZ 


s D</z 


(I+«,^*)(|+2U.2) 

The  left-hand  side  of  (194)  can  be  separated  into  three  terms,  and 
becomes: 


'i2dZ  ^ zu-dZ  2U-JZ 

By  integration: 


ifzu2 


(191) 


(192) 

(193) 

(194) 
(194) 

(195) 


(196) 


F 


f 


"mmmm 
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where  G is  the  integration  constant.  Eq.  (196)  can  be  rewritten  as: 

zPdWjz 

-Ufluicto/w  2Z 

K 

substituting  (191)  in  (197)  and  considering  that  D = 


(197) 


2Bb(l+u‘) 


2 — , we  find: 


p o . nil  r-  R. 


(198) 


With  the  initial  conditions  Pj(0)  = P^,  Pjj(O)  = 0,  we  find  that  G = 1 and, 
finally: 


p..p.+  iu;5^  . 


(199) 


Eq.  (199)  can  be  rewritten  as: 


P - ? 2 

'IT  ~ 


R JJ 


(200) 


(201) 


4.5  Discussion  of  Solution. 

Eqs.  (200)  and  (201)  show  that  P,^  increases  toward  a maximum  of 
P P 

(1  + -===),  and  P_  decreases  toward  a minmam  of  (1  — ) , where  P^^ 
/l+u2  ^ ^ <^l+u2  “ 

reaches  its  maximum.  Pj  reaches  its  minimum  when  (z  -►  ») : this  incorrect 

result  is  a consequence  of  the  fact  that  (175)  and  (176)  correspond  to 


increasing  Pjj  and  decreasing  Pj. 
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According  to  relations  (181)  and  (182),  we  must  switch  to  equations  (178) 

2 P 

and  (179)  when  P..  reaches  the  value  P^  this  value  is  smaller  than  -^(1+-===). 

“ ‘^1+u'^  dPjj  M+u2 

Also,  we  know  that  when  P^t  increases,  — decreases;  one  bounce 

’’o  ''o“^ 

before  the  switch  takes  place  we  have  that  P.  - =■  and  P^^  - =r,  the  values 

dPjj  1+u  1+u 

of  Pj  and  Pjj  for  which  ■ is  minimum.  Substituting  in  (176)  we  find: 


R,tf 


+ 2UL 


P 


(202) 


The  correct  solution  is  therefore  obtained  in  the  following  manner.  P,(z) 
and  Pjj(z)  must  follow  formulas  (200)  and  (201)  until  ^ per  bounce 

length;  for  the  next  bounce,  we  must  switch  to  equations  (178).  and  (179),  whose 
solutions  are  obtained  from  (200)  and  (201)  by  interchanging  P,  and  P.,,.  IVhen 

V"'  '■o 

the  switch  occurs,  P,,  is  greater  than  j smaller  than  j. 

1+u  1+u 

Our  solution  is  shown  in  Fig.  6,  where  the  ratio  PjjC2)/Pq  is  plotted  as 
a function  of  z for  different  values  of  the  parameter  u.  In  the  same  figure 
we  also  plotted  the  Miller-Marcatili  solution,  as  given  by: 

Pjj  = Pg  cos^  Cz  (203) 


I 


with 


SVv~t,C 

?pv>  u-u>  2Kj'<s 


(204) 


We  see  that  the  Miller-Marcatili  solution  approaches  ours  as  u increases,  i.e. 
as  c +■  <*>.  In  Fig.  6 the  curves  have  been  stopped  before  Pjj  begins  to  decrease. 

According  to  our  solution,  there  cannot  be  complete  power  transfer  from 
guide  I to  guide  II.  The  maximum  power  transfer  possible  increases  with  the 


o 
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Tf 

Figure  6.  Normalized  power  in  guide  II  as  a function  of  the  coupling  length. 

Our  exact  solution.  Marcatili's  approximation. 


F 


separation  distance  c,  and  it  reaches  100  percent  in  the  limit  c = 

For  c = 0,  the  two  guides  form  one  guide  with  width  2b,  then  u = 0 and 
from  (176)  we  find  that 

Pn  PoC'-ePb  ) ^ (205) 

which  means  that  after  a few  bounce  lengths  from  z = 0,  almost  half  the 

incident  power  is  transferred  to  guide  II.  The  maximum  transfer  of  power 
p 

to  guide  II  is  In  fact,  the  powers  in  guides  I and  II  should  be  equal 

after  the  first  bounce.  We  do  not  obtain  this  result  because  our  solution 
is  valid  only  for  a small  amount  of  power  transfer  per  bounce. 
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CONCLUSION 

In  our  procedure  for  solving  the  coupling  between  dielectric  optical 
waveguides  we  have  used  three  steps: 

STEP  1:  Calculate  the  reflection  and  transmission  coefficient  at  the  wave- 

guide boundaries.  We  saw  that  phase  and  magnitude  of  these  two  coefficients 
are  functions  of  the  separation  distance  between  the  guides. 

STEP  2:  Solve  the  eigenvalue  equation  for  modes  in  the  region  (0  - z - L) . 

We  showed  that  for  not  too  small  values  of  the  separation  distance,  th'  inci- 
dent mode  in  guide  I in  the  region  (z  - 0)  couples  into  a single  mode  in 
the  region  (0  - z - L) , and  the  two  modes  have  the  same  propagation  constant. 

With  this  condition,  we  may  neglect  reflections  at  the  boundaries  (z  = 0) 
and  (z  s L) . 

STEP  3:  Derive  the  differential  equations  describing  the  power  transfer  between 

the  two  guides  in  the  most  general  case.  All  the  coefficients  in  these  equa- 
tions have  been  explicitly  calculated.  Because  of  the  fact  that  the  trans- 
mission coefficient  between  the  two  guides  leads  the  reflection  coefficient 
by  the  phase  difference  between  the  fields  in  the  two  waveguides  is  a 
constant  along  z,  and  equals  - With  this  condition,  the  solution  of  the 

differential  equations  for  the  power  transfer  has  been  easily  obtained.  How- 
ever, the  form  of  the  solution  is  not  simple.  We  saw  that  there  are  two 
different  solutions,  one  valid  when  the  power  Pjj  in  guide  II  is  increasing, 
the  second  when  Pjj  is  decreasing.  However,  rules  describing  how  to  switch 
from  the  first  solution  to  the  second  were  found;  this  switch  occurs  when  P^j 
reaches  its  maximum.  We  showed  that  there  is  complete  transfer  of  power  only 
when  the  separation  distance  tends  to  infinity  (extremely  weak  coupling) . 

Although  our  theory  was  developed  for  the  particular  case  of  rectangular 
waveguides,  it  could  be  extended  to  waveguides  of  different  cross  sections.  For 
example,  it  could  even  be  used  to  describe  the  coupling  of  two  optical  fibers  via  a 
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coaxial  dielectric  sleeve;  this  specific  problem  will  be  studied  in  a 
forthcoming  paper. 

Finally,  we  observe  that  our  theory  allows  us  to  use  the  directional 
coupler  as  a mode-selecting  device.  If  several  modes  are  incident  in  guide  I, 
we  can  choose  the  separation  distance  c and  the  coupling  length  L to  maximize 
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